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Abstract. In this paper we study geometric, algebraic, and computational aspects of 
flexibility and infinitesimal flexibility of Kokotsakis meshes. A Kokotsakis mesh is a mesh 
that consists of a face in the middle and a certain band of faces attached to the middle 
face by its perimeter. In particular any (3 x 3)-mesh made of quadrangles is a Kokotsakis 
mesh. We express the infinitesimal flexibility condition in terms of Ceva and Menelaus 
theorems. Further we study semi- algebraic properties of the set of flexible meshes and 
give equations describing it. For (3 x 3)-meshes we obtain flexibility conditions in terms 
of face angles. 
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Introduction 

Consider a mesh having one special n-vertex face and a "strip" of faces attached to 
the special face along the perimeter, such that any vertex of a special face is a vertex of 
exactly three faces in the strip. This mesh is called a Kokotsakis mesh. Kokotsakis meshes 
with polygonal faces are known in the literature as polyhedral surfaces of Kokotsakis type. 
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In this paper we study conditions of flexibility and infinitesimal flexibility for Kokotsakis 
meshes. 

The first steps in studying such meshes were made in f 932 by A. Kokotsakis in his pa- 
per [5], where he found conditions of infinitesimal flexibility. Recently Kokotsakis meshes 
for 4-gons (such meshes are also called (3x3) -meshes) were used in the study of discrete 
conjugate nets flexibility. A. I. Bobenko, T. Hoffmann, W. K. Schief proved in [3] that 
a nondegenerate discrete conjugate net is isometrically deformable if and only if all its 
3x3 subcomplexes are isometrically deformable (3x3)-meshes. In [3] they also showed 
that the conjugate nets admit infinitesimal second-order deformations if and only if their 
reciprocal-parallel surfaces constitute discrete Bianchi surfaces (see also [2]). 

In this paper we give a geometrical interpretation of infinitesimal flexibility in terms 
of Ceva and Menelaus configurations. Further we work with equations that define the 
set of flexible Kokotsakis meshes, in particular we show semi-algebraicity of this set. 
The technique developed in this paper allows us to calculate linear independence of the 
gradients of the first 6 flexibility conditions for a chosen flexible quadrangular Kokotsakis 
mesh. This means that we can check if the codimension of a stratum of meshes is smaller 
than or equal to 6. The computational complexity of the problem does not allow us to 
proceed further. We remind the reader that the codimension of the set of flexible Voss 
meshes as well as that of symmetric flexible meshes is 8, for information on Voss surfaces 
we refer to [3]. So it is most probable that the codimensions of irreducible components 
coincide, and equal 6, 7, or 8. If we are lucky and these codimensions are exactly 8, 
than the geometric solution of the problem can be obtained by enumeration of all possible 
components. 

To check the conditions at a single point we introduce more simple "flexibility vector 
calculus" conditions. In the case of (3x3)-meshes we now able to calculate the 8 first such 
conditions. 

Finally, for the case of quadrangular Kokotsakis (3x3)-meshes we show that the prop- 
erty of the mesh to be flexible depends only on the shapes of faces but not on the geometry 
of the mesh in IR 3 . We propose an algorithm to write these formulae explicitly. 

Organization of the paper. We start in Section 1 with showing an infinitesimal 
flexibility condition in algebraic and geometric terms. Section 2 studies properties of the 
set of flexible meshes. We describe equations defining this set. In Section 3 we introduce 
a technique to write flexibility conditions in terms of angles of the faces of meshes. We 
conclude the paper in Section 4 with some open questions for further study. 

General notation. To avoid the annoying description of multiple cases we consider 
the index i belonging to the cyclic group Z/nZ. We also choose a partial ordering 1 < 
2 < 3 < . . . < n for this group. 

We denote the vertices of the special face in the middle of the mesh by i — 1, . . . , n, 
preserving their order at the boundary. Let Oj = A i+1 — A^. Notice that by definition we 
have J27=i °i = 0- 
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Figure 1. Notation. 



Consider a vertex Ai (i — 1, . . . , n). We denote vertices V* and Wi such that the four 
faces meeting at Ai contain angles A i+ iAiAi-i, A^iAiVi, ViAiWi, and WiAiA i+i . We 
denote these angles by at, fa, 7$, and </?j, respectively. Let t> j = — and iUj = Wi — A^. 

Remark. Without loss of generality we assume that the special face A\A 2 . . . A n is fixed, 
while the other 2n faces may move. 

Denote by (ui,u 2 ) and by [1*1,1x2] the scalar and vector products, respectively. 



1.1. Algebraic conditions. In this subsection we recall some results on infinitesimal 
flexibility from the work [5] by A. Kokotsakis. We rewrite his results in a slightly modified 
way more convenient for us. 

Consider a (2x2)-mesh with a common vertex A. Let us fix a face BAC. Let BAV, 
VAW, and WACbe the other three faces. Denote by b, c, v, and w the vectors B — A, 
C — A, V — A, and W — A respectively. This mesh is flexible in general position and has 
one degree of freedom. 

Proposition 1.1. Suppose that the vectors v and w are not collinear and the vectors b 
and c are not in the plane spanned by the vectors v and w. Then for any infinitesimal 
motion (v,w) of the described (2x2) -mesh there exists a non-zero real A such that 



Proof. Any infinitesimal motion is a solution of the following linear system with 6 variables 
(coordinates of v and w) and 5 equations: 



1. Infinitesimal flexibility of Kokotsakis meshes 




(v,v) = 



(w, w) — 



< (v, w) + {v, w) — . 



(b,v) = 
k (c,w) = 
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On one hand the solution mentioned in the statement of the proposition satisfies this 
system. On the other hand, if the conditions on vectors b, c, v, and w hold then the rank 
of the linear system is 1 (this follows from direct calculations in coordinates). □ 

From Proposition 1.1 we immediately have the following statement (this statement is 
another formulation of the result in [5] by A. Kokotsakis). For a Kokotsakis mesh M put 
by definition 

(a,i-i,Vi,Wi) 



x(m) = n 



| K Oi,Vi,Wi) 

Corollary 1.2. A Kokotsakis mesh M is infinitesimally flexible if and only ifx(M)—l. 

□ 

1.2. Geometric infinitesimal condition for a Kokotsakis mesh with planar faces. 

Consider points Ai, . . . , A n in a plane. We denote this plane by ir. Let U denote the 
intersection of the plane containing the points Ai, Vi, and Wi with the plane n. We 
suppose that the lines U and are not parallel, and denote their intersection point by 
B t . 

For a (3x3)-mesh constructed with planar faces the infinitesimal flexibility condi- 
tions are shown in the following proposition. For more information see the work [6] 
by H. Pottmann, J. Wallner and the book [7] by R. Sauer. 

Proposition 1.3. Consider a mesh in general position. Then it is infinitesimally flexible 
if and only if the lines B 1 B 3 , A 2 , A 3 , and A 4: A 1 intersect in one point (i.e. the points A i} 
B i7 i — 1,2,3, 4 are aligned in a Desargues configuration). 

Let us generalize this statement for an arbitrary Kokotsakis mesh with planar polygons 
as faces. For i — 1, . . . , n we define real numbers t; L from the following afline equations: 

Ai = tiBi_ x + (1 - ti)Bi. 

Theorem 1.4. The closed broken line AiBiA 2 B 2 . . . A n B n A\ satisfies the generalized 
Ceva-Menelaus condition, i.e. 

n 

nV 1 =(-ir- 

Remark 1.5. If there is one point Bi "at infinity" one should replace 

1 ~ t i . 1 ~ ^+1 ^ ( 1) £ ~ ti+1 

ti ti 

in the product. Here ti and 1— t i+ i denotes the absolute values of the vectors Ai—B^ and 
B i+ i—A i+1 respectively. We take e = if the directions of these two vectors are opposite 
and e — 1 if they coincide. 

We leave for the reader the other cases of some of the vertices B, L being "at infinity" as 
an exercise. 
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Proof. Choose some orientation on the plane ir. For an ordered couple ei,e 2 of vectors 
denote 

{1 if (ei, e-i) is a positively oriented basis of the plane 
if e\ and e 2 are collinear 
—1 otherwise 

Denote by hpQR the (non-oriented) distance from the point P to the line QR. Then 
for any i G {1, . . . , n} we have: 

(vi, Wi, Oi_i) sgn( J B i _i5j, 
sgn(S i _iS i ,i4iA+i) |AA+i|/sin(A+iA^) 



Therefore, 



sgn(S i _iS i ,i4i_iA) |A-iA|/sin(A-iAP*-i)' 

n 1 ~^ _ TT sgn(j?i-2-Bi-i, A-iA)|Pj-iA| 
:1 *< "fj -sgn(S iJ B m ,AA+i)|A^| 
sgn(Sj_iSj, M+i) |AA+i| sin(A-iA-B 



n 



n 



i-ly 



* \sgn(Sj_iSj, A-iA) |A-iA| sm(A+iA-B* 



C_1"| n TT \ v ii w ii a i) /_i\n 
1 J '11 Oi-i) 1 J • 

The second equality holds by the law of sines for the triangles AiA i+ iBi for i — 1, . . . , n. 

□ 

Denote by [A, B; C, D] the point of intersection of the lines passing through A, B and 
C, D respectively. 

Denote P = B 1 . Let Oi = [B n , A i+1 ; B i+1 , P;_i], and Pi = [B i+1 , B n ; B h Oi\, for 
i = 1, . . . , n— 2. 

Proposition 1.6. Geometric condition I. The condition 

i=i 4 

?5 equivalent to the condition that the points A n and P n -2 coincide {see Figure 2.1). 

In the case n = 3 we have: £/ie /ines -BiA; B 2 A 1; and B 3 A 2 intersect at one point (Ceva 
configuration). In the case n = 4 we have: the lines BxB 3 , A 2 ,A 3 , and A 4: A 1 intersect at 
one point (Desargues configuration) . 

Proof. We prove the proposition by induction on n. 

By Ceva's theorem and Theorem 1.4 the statement holds for all configurations with 
n — 3. For pairwise distinct collinear points A, B, C denote 

i a r>r*\ — / 1 ^ ^ ne P°i n t ^ i s m the segment [A, C] 
^ ' ^ — 1 otherwise 
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I: The points ^5 and P3 coincide. 



II: The points A4, A5, and P2 are in a line. 



Figure 2. Geometric conditions I and II for a pentagon. 



Suppose the statement holds for all configurations with n — 4. By Ceva's Theorem in 
the triangle B n B 1 B 2 we have: 



AB t 



i=l ' ' x ' i=3 

Then by induction the statement is true for the broken line P1B2A3B3 . . . A n B n P\. □ 
Proposition 1.7. Geometric condition II. The condition 



n 

i=i 



1 -t: 



-1 



zs equivalent to the condition that the points P n s, A n _\, and A n are collinear {see Fig- 
ure 2. IT). 

In the case n = 3 we have: the points A x , A 2 , and A 3 are contained in one line (Menelaus 
configuration). 

Proof. We prove the proposition by induction on n. By Menelaus' theorem and Theo- 
rem 1.4 the statement holds for all configurations with n — 3. The induction step is the 
same as in the proof of Proposition 1.7, so we omit it here. □ 

Remark 1.8. The condition of Corollary 1.2 can be rewritten as Geometric condition I in 
the case of even n, and as Geometric condition II in the case of odd n. 

Corollary 1.9. Any Kokotsakis mesh with a triangle in the middle, in general position, 
is not infinitesimal flexible (and therefore it is not flexible) . □ 

In particular we have the following conclusion. Consider a closed mesh in general 
position (i.e. there are no 3 vertices lying in one line) whose edge-graph is 4-valent. Then 
the mesh is not flexible if at least one of its faces is a triangle. Actually this means that the 
Euler characteristic of the surface is non-positive. So a flexible mesh in general position 
is not homeomorphic to the sphere or to the projective plane. 
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2. Flexibility conditions 



2.1. Semialgebraicity of the set of flexible Kokotsakis meshes. First we define a 
configuration space of all Kokotsakis meshes whose central face is an n-gon. Any Kokotsa- 
kis mesh contains 3n points. So the configuration space of all possible vertex configurations 
is (R 3 ) 3r \ We will not consider configurations for which some face has three consequent 
vertices in a line, or which has faces with coinciding vertices. Denote this set by E n . Fi- 
nally, the configuration space of all Kokotsakis meshes is by definition (R 3 ) 3n \ E n , denote 
it by n . 

Notice that we are not restricted to the case of planar faces in this section. Consider 
the vector field: 



and denote it by £„. This vector field is clearly everywhere defined in n . 

Further when, in this section, we speak of a motion we assume that the special face in 
the middle is fixed while the others are moving. 

Theorem 2.1. Consider a flexible mesh in n . The motion trajectory of the mesh in 
some small neighborhood is contained in some integral curve of the vector field £ n . 

Proof. By Proposition 1.1 a tangent vector for the trajectory of a flexible mesh M in O n 
is proportional to the vector £(M). This implies the statement of the theorem. □ 

Denote the set of all flexible meshes by -Fe„. Let us consider an infinite sequence fi of 
algebraic functions on (IR 3 ) 3ra , defined inductively: 
%). Take f\ to be the numerator of the rational function x(M) — 1. 
ii). Take f n to be the numerator of the rational function (grad(/ n _i), £ n ). 

Remark 2.2. We ignore denominators since they are nonzero algebraic functions on n . 

Theorem 2.3. The closure of F® n in 6 n is a semialgebraic set coinciding with 



Proof. Notice that the vector field £(M) is non-zero at all the points of © n . 

By Proposition 1.1, for a flexible mesh M not contained in E n the velocity vector is 
proportional to the vector £(M), and therefore satisfies all the conditions f n = 0. 

On the other hand all functions f n are algebraic and therefore their common locus is an 
algebraic variety (by Hilbert's basis theorem it is sufficient to take only a finite number of 
such equations). Therefore in a nonsingular point of this variety we have all the existence 
conditions for a finite motion along the field £(M) (see in [1]). Since the set of singularities 




for i — 1, . . . , n. 



and 



i=i 
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of an algebraic variety is nowhere dense in this variety, the closure of the set of nonsingular 
points of a variety is the variety itself. □ 

2.2. Vector calculus of flexibility conditions. In practice it is quite hard to calculate 
the above /„ explicitly. To check some given configuration we recommend a slightly 
different technique. 

Since x(M) = along any trajectory, we immediately have many conditions of flexi- 
bility: x^ n \M) = 0, where the derivatives are taken with respect of the time parameter 
of the vector field. We recommend to simplify the obtained function after each round of 
differentiation using previous ones. 

Let us introduce basic ingredients to calculate derivatives. We use a prime to indicate 
differentiation. 

(a i:Vi y = -( fli _ 1 , Ql , ^n^gf; 

i=i 

i-l 



(a^^a^Vi)' = ({ai^a^^Vi) - (a^i, a i )(a i _i, J[ (a ( J a .^ w . . 
(ai-i,ai,Wi)' = ({ai-i,ai){ai,Wi) - (oj, a i )(a i _i,«; i )) fl 
(a fc ,^,Wi)'= ((a fc ,Ui)(ai_i,«;i) - (ai_i,a fc > J fl %S 
{a h a k ){vi,Wi) - (ai,Vi){a k ,Wi)) f\ 



i-l 

i , . f j j i i 



In fact, to calculate x (n) ( M ), 

one only needs the cases k — i— 1 and A; = i for the last 

formula. 

Example 2.4. The second condition (%(M))' = is equivalent to 

/ (ai,^i)(a,-i,w,) - (a,-i,ai)(i;i,Wt) irr (flj-i,^, Wj) j 

(ai,ai)(vi,Wi) - (a h Vi)(a h Wi) yj (aj-^v^Wj) 



\ ^ "i/ y^j, Wj/ - Ujl \ U U W H TT 

^ \ (oi,Vi,Wi) Al (aj,Vj,Wj) 

A / (ai-i,^)(ai_i,Wi) - (a^-i, a^_i)(^, w^) 4-r (a 3 -_i, Wj) 

| (Qi-i,ai)(^,Wi) - (a i ,t; i )(a i _i,Wi) -pi- (a^i, ^ \ _ 
^ I (ai_i,^,^) Al {aj,Vj,w-) J 
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3. Internal flexibility conditions for (3x3)-meshes 

In this section we describe a technique to verify the flexibility of a (3x3)-mesh with 
planar faces in terms of inner geometry of a mesh, i.e., in terms of the fixed angles between 
the edges of the same faces. 

Consider a (3x3)-mesh with a planar quadrangle in the middle. Denote the 

angles between the faces adjacent to AiA i+ i by u>i (i — 1, . . . , n). 

Lemma 3.1. The following relation holds for any vertex A { 

cos «j cos fa cos <fi + sin aj sin fa cos ipi cos Ui_i + sin ccj cos fa sin ipi cos u— 
cos aj sin fa sin ipi coscjj„i cosc^j + sin fa sin <pi sina;j_i sinwj = cos7j. 

Denote this equation by (Ri). 

Proof. Note that the values of the angles do not depend on the lengths of vectors a i: Oj_i, 
Vi, Wi. So we can choose them equal 1. Let us make calculations in the coordinate system 
where = (1,0,0), and Oj_i = (cos a, sin a, 0). Then we have 

Vi = (cos (pi, sin ipi coscjj, siny? sincjj), 

= cos/3j(coso;j, sinctj, 0) + sin/3j ^ coscjj_i(sin a iy — cosaj, 0) + sinwj_i(0, 0, l)j 

The formula of the lemma coincides with (vi,Wi) = cos7j and therefore holds. □ 
We get the desired equations in 4 steps. 

Step 1. Denote cosc^j by t; L . The first step is to get rid of sincUj sinu;j_i in formula (Ri). 
This can be done by putting the corresponding term to the right and the rest to the left, 
than squaring both sides, and finally substituting sin 2 oji sin 2 cjj_i = (1 — tf)(l — t 2 ^). 
Denote the resulting equation by (Ri). It leads to 

Cl,i + C^iti-i + Cz,iti + C , 4 j jt 2 _ 1 + C^iti-]ti + C*6,itf + C , 7 i jt 2 _ 1 tj + Cs,iti-it1 + Cg i jt 2 _ 1 t 2 = 0, 

where 

Cl,i — ( cos a i cos fa cos V 9 * — cos li) 2 ~ s i n2 A srn2 ¥u 

C 2; i = 2(cosaj cos fa cos(/9j — cos7j) sinctj sin fa cos ipf 
C 3;i = 2(cosaj cos fa cos ipi — cos7j) sin«j cos fa sin ipf 
C^i = sin 2 ai sin 2 $ cos 2 ipi + sin 2 sin 2 ipf 

C^i = 2(cos fa cos (fi — 2 cos 2 cos $ cos </?j + cos ai cos 7$) sin fa sin 
C 6j j = sin 2 ai cos 2 /3j sin 2 ipi + sin 2 sin 2 iff 
Cj^ = —2 sin ai cos ctj sin 2 sin tpi cos 
6*8,1 = —2 sin otj cos «j sin $ cos $ sin 2 
69,1 = — sin 2 ai sin 2 $ sin 2 ipi. 

Step 2. Now we recursively evaluate tj as a linear function in U whose coefficients are 
rational functions of U-\. We will need only the formulae for t\ and t\. 

Step 3. We use (R3) and the formula for t\ given in Step 2 to evaluate t 2 as a function 
of ti and £ 3 , and analogously we employ (R4) with the formula for t\ given in Step 2 to 
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evaluate t 4 as a function of t 1 and t 3 . We get 

. (C , 9,2^+C8,2*l+C6,2)(C 6 ,3ti+C'3,3t3+Cl,3)-(C9,3*i+C7,3t3+C4,3)(C'4,2tf+C2,2*l+Cl,2) . 

2 (C , 9,2t|+C , 8 ,2il+C 6 ,2)(C 8 ,3t|+C 5 ,3t3+C2,3)-(C , 9,3t|+C 7 ,3i3+C4,3)(C' 7 ,2^+C5,2tl+C3,2)' 

< (C9,l +C7,if i +0*4,1 ) (£4,4^3+02, 4*3 +Cl, 4) — ((^9,4*3 +Cs, 4<3+C6, 4)(C , 6,lt- L +C3,l*l+Cl,l) 

4 ~ (C9,l^^+C7,l^l+C4,l)(C 7 ,4^+C5,4^3+C3,4)-(C*9,4^+C* 8 , 4 ^3+C6,4)(C*8,l^^+C , 5,l ^ l+C'2,l)• 

.Step 4- Now for ti and t 3 we have the system 

f Ri(ti,U(t 1: t 3 )) = 
1 R2(t 1 ,t 2 (t 1 ,t 3 )) = ' 

First, note that the functions i?i and i?2 are rational in variables t\ and £ 3 . So we can 
restrict ourselves to consider the numerators. Denote these numerators by fi(ti,t 3 ), and 
f2(ti,h)- The functions f\ and / 2 are polynomials of degree 10, and we need to check if 
these two polynomials define a non-discrete set. This happens only if f\ and f'2 have a 
common nonconstant factor. Take a resultant of f\ and f'2 considered as polynomials in 
t 3 . This resultant denote it by p{t\) is a polynomial in t\. 

Theorem 3.2. A mesh in ©„ is flexible if and only if all the coefficients of the polynomial 
p(ti) equal zero. 

Proof. We give an outline of the proof. 

It is clear that if a mesh is flexible then the set of common solutions of Ri = for 
i = 1,2,3,4 is nondiscrete and therefore one of the equations can be eliminated. By 
Proposition 1.1 in the case of meshes in n all the t £ 's are nonconstant while the mesh is 
moving, so the resultant p(£i) is identically zero. 

Suppose now that all the coefficients of the polynomial p(ti) equal zero, then there are 
at most three independent equations among R{ = for i = 1,2,3,4. The common set of 
solutions of the equations Ri = is almost always at least one-dimensional. 

The dimension of this set is zero only if the corresponding mesh is a limit of a sequence 
of flexible meshes with possible motion curves shrinking to this mesh. Such mesh is 
infinitesimally flexible with more than one degrees of freedom. By Proposition 1.1 such 
meshes are not in 0„. □ 

Corollary 3.3. Flexibility of meshes in 9 n is determined by the shapes of faces (but not 
of an embedding in the space as a mesh) . □ 

To simplify symbolic calculations we recommend to compute the resultants of f\(n,t 3 ) 
and f2{n,t 3 ) for —50 < n < 50. If all 101 (=10 2 + 1) resultants evaluate to zero, than the 
polynomial p{ti), which has at most 100 zeroes, has all zero coefficients. 

Finally, we give an algorithm that calculates the functions fi(x,y) and f 2 (x,y) (as f[l] 
and f[2] in the output below). 

TwoEquations:=proc (alpha, beta, gama, phi, N) local i,C,R,t,f ,tn,td; 

for i from 1 to 4 do 

C [1 , i] : = (cos (alpha [i] ) *cos (beta [i] ) *cos (phi [i] ) -cos (gama [i] ) ) "2- 
sin(beta[i] )~2*sin(phi [i] )~2; 
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# Here we write all the formulae for C[2,i] , . . . ,C[9,i] from Step 1 
end do : 

tn [2] : =- ( (C [6 , 2] +C [8 , 2] *t [1] +C [9 , 2] *t [1] ~2) * (C [1 , 3] +C [6 , 3] *t [3] ~2+C [3 , 3] *t [3] ) - 
(C [7 , 3] *t [3] +C [4 , 3] +C [9 , 3] *t [3] ~2) * (C [1 , 2] +C [2 , 2] *t [1] +C [4 , 2] *t [1] ~2) ) : 
td [2] : = ( (C [6 , 2] +C [8 , 2] *t [1] +C [9 , 2] *t [1] ~2) * (C [8 , 3] *t [3] ~2+C [5 , 3] *t [3] +C [2 , 3] ) - 
(C [7 , 3] *t [3] +C [4 , 3] +C [9 , 3] *t [3] ~2) * (C [3 , 2] +C [5 , 2] *t [1] +C [7, 2] *t [1] "2) ) : 

tn [4] : =- ( (C [1 , 4] +C [2 , 4] *t [3] +C [4 , 4] *t [3] ~2) * (C [4 , 1] +C [7 , 1] *t [1] +C [9 , 1] *t [1] "2) - 
(C [6 , 4] +C [8 , 4] *t [3] +C [9 , 4] *t [3] ~2) * (C [1 , 1] +C [3 , 1] *t [1] +C [6 , 1] *t [1] "2) ) : 
td [4] : = ( (C [3 , 4] +C [5 , 4] *t [3] +C [7 , 4] *t [3] ~2) * (C [4 , 1] +C [7 , 1] *t [1] +C [9 , 1] *t [1] '2) - 
(C [6 , 4] +C [8 , 4] *t [3] +C [9,4] *t [3] ~2) * (C [2, 1] +C [5 , 1] *t [1] +C [8, 1] *t [1] ~2) ) : 

RBar [1] : = (C [1 , 1] +C [6 , 1] *t [1] ~2+C [3 , 1] *t [1] ) *td [4] ~2+ 

(C [8 , 1] *t [1] ~2+C [5 , 1] *t [1] +C [2 , 1] ) *td [4] *tn [4] + 

(C [7 , 1] *t [1] +C [4 , 1] +C [9 , 1] *t [1] ~2) *tn [4] ~2 : 
RBar [2] : = (C [1 , 2] +C [2 , 2] *t [1] +C [4 , 2] *t [1] ~2) *td [2] "2+ 

(C [3 , 2] +C [5 , 2] *t [1] +C [7 , 2] *t [1] ~2) *td [2] *tn [2] + 

(C [6 , 2] +C [8 , 2] *t [1] +C [9 , 2] *t [1] ~2) *tn [2] ~2 : 

f [1] : =(x,y)-> (subs ({t [1] =x,t [3] =y}, simplify (RBar [1]))) : 
f [2] : = (x , y) -> (subs ({t [1] =x , t [3] =y} , simplify (RBar [2] ) ) ) : 

return (f ) : 
end proc; 

Now we can compute resultants of f[l](yo,y) and f[2](y ,y) for any value y with any 
precision N, or symbolically. 

yO:=l; 
N:=50; 

f : =TwoEquations (alpha, beta, gama, phi ,N) : 
ff [1] :=evalf [N] (f [1] (yO,y)) : 
ff [2] :=evalf [N] (f [2] (yO,y)) : 

evalf [N] (resultant (f f [1] /coeftayKf f [1] ,y=0,4) ,f f [2] /coeftayKf f [2] ,y=0,4) ,y)) ; 

Still we do not know which of the listed 101 equations are independent, the first inter- 
esting problem in this direction is the following. 

4. Some questions for further study 

4.1. Geometry of conditions. In Section 1 we studied the geometric conditions for 
infinitesimal rigidity of meshes with planar faces, the following question arise here. 

Problem 1. What is the geometric infinitesimal flexibility condition in case of meshes 
with nonplanar faces? 

A similar problem remains open for the case of flexibility conditions introduced in 
Section 2. 

Problem 2. Find a geometric interpretation of the condition y/ n )(M) = for n > 1. 
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We do not know how the condition 'x'(M) = is related to the second-order rigidity 
condition defined by A. I. Bobenko, et al. in [3]. 

4.2. Codimension calculation. In Section 2 we described the closure of all flexible 
meshes F® n by algebraic equations. We conjecture that it is a complete intersection. It 
is most probable that all the components of the closure are of the same dimension. 

Remark 4.1. Let us say a few words about the case of (3x3)-meshes. Direct calculation 
of the tangent space in the points of Fq 4 based on the results of the second section shows 
that the codimension of some of the irreducible components of F® A is greater than or 
equal to 6. On the other hand, the Voss surfaces (c.f. [8] and [3]) forms a set of flexible 
meshes of codimension 8. In general we suspect that all the irreducible components of 
Fq a has the same codimension between 6 and 8. Note that for codimension 8 we know 
the following families of flexible meshes: 

i) Voss surfaces (defined by equations ctj = 7,, fa — <fi for i — 1, 2, 3, 4). 

ii) symmetric meshes where the angles at the vertices A 1 and A 2 coincide with the 
angles at the vertices A 4 and A 3 : 

"i = "4, Pi = V?4, 7i = 74, (Pi = Pa] oc 2 = a 3 , p 2 = <p 3 , 72 = 73, ^2 = th- 
in) symmetric meshes where the angles at the vertices A 1 and A A are the same as at 
the vertices A 2 and A 3 . 

iv) meshes obtained from the cases i) — iii) by replacing some Vi and/or Wj by —Vi 
and/or —Wj. For instance, if we replace v\ by —v\ in the case of Voss meshes, then we 
get equations ot\ + 71 = n, Pi + (fii — ir, and a« = 7$, Pi = ipi for i — 2,3, 4. 

If the codimension of F@ A is 8, then the natural problem is to complete the list of flexible 
meshes in general position. 

All this leads to the following problem. 

Problem 3. Find the codimension of the closure of F e „ for n — 4, 5, . . . 

The same question is interesting for the case of an algebraic set defined by 101 equations 
in Section 3, does this set has the same codimension? 
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